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ABSTRACT 


The  deviatoric  and  isotropic  source  components  for  17  explosions  at  the  Nevada  Test  Site,  as  well  as  12  earthquakes 
and  3  collapses  in  the  surrounding  region  of  the  western  US,  are  calculated  using  a  regional  time-domain  full 
waveform  inversion  for  the  complete  moment  tensor.  The  events  separate  into  specific  populations  according  to 
their  deviation  from  a  pure  double-couple  and  ratio  of  isotropic  to  deviatoric  energy.  The  separation  allows  for 
anomalous  event  identification  and  discrimination  between  explosions,  earthquakes,  and  collapses.  Confidence 
regions  of  the  model  parameters  are  estimated  from  the  data  misfit  by  assuming  normally  distributed  parameter 
values.  We  investigate  the  sensitivity  of  the  resolved  parameters  of  an  explosion  to  imperfect  Earth  models, 
inaccurate  event  depths,  and  data  with  a  low  signal-to-noise  ratio  (SNR)  assuming  a  reasonable  azimuthal 
distribution  of  stations.  In  the  band  of  interest  (0.02-0.10  Hz)  the  source-type  calculated  from  complete  moment 
tensor  inversion  is  insensitive  to  velocity  model  perturbations  that  cause  less  than  a  half-cycle  shift  (<5  sec)  in 
arrival  time  error  if  shifting  of  the  waveforms  is  allowed.  The  explosion  source-type  is  insensitive  to  an  incorrect 
depth  assumption  (for  a  true  depth  of  1  km),  but  the  goodness-of-fit  of  the  inversion  result  cannot  be  used  to  resolve 
the  true  depth  of  the  explosion.  Noise  degrades  the  explosive  character  of  the  result,  and  a  good  fit  and  accurate 
result  are  obtained  when  the  SNR  is  greater  than  5.  We  assess  the  depth  and  frequency  dependence  upon  the 
resolved  explosive  moment.  As  the  depth  decreases  from  1  km  to  200  m,  the  isotropic  moment  is  no  longer 
accurately  resolved  and  is  in  error  between  50%  and  200%.  However,  even  at  the  most  shallow  depth  the  resultant 
moment  tensor  is  dominated  by  the  explosive  component  when  the  data  has  a  good  SNR.  The  sensitivity 
investigation  is  extended  via  the  introduction  of  the  network  sensitivity  solution,  which  takes  into  account  the  unique 
station  distribution,  frequency  band,  and  SNR  of  a  given  test  scenario.  An  example  of  this  analysis  is  presented  for 
the  North  Korea  test,  which  shows  that  in  order  to  constrain  the  explosive  component  one  needs  a  certain  station 
configuration.  In  the  future  we  will  analyze  the  bias  in  the  source-type  parameters  due  to  error  in  the  Green’s 
function  by  incorporating  a  suite  of  suitable  velocity  models  in  the  inversion. 
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OBJECTIVES 

We  implement  the  time-domain  full-waveform  inversion  of  regional  data  for  the  complete  moment  tensor  devised 
by  Minson  and  Dreger  (2008)  after  Herrmann  and  Hutcheson  (1993)  based  on  the  work  of  Langston  (1981).  In 
general,  synthetic  seismograms  are  represented  as  the  linear  combination  of  fundamental  Green's  functions  where 
the  weights  on  these  Green's  functions  are  the  individual  moment  tensor  elements.  The  Green's  functions  for  a  one¬ 
dimensional  (1-D)  velocity  model  of  eastern  California  and  western  Nevada  (Song  et  al.,  1996)  are  calculated  as 
synthetic  displacement  seismograms  using  a  frequency-wavenumber  integration  method  (Saikia,  1994).  The 
synthetic  data  are  filtered  with  a  4-pole  acausal  Butterworth  filter  with  a  low  comer  of  0.02  Hz  and  a  high  comer  of 
0.05  Hz  and  0.1  Hz  for  events  with  Mw  >  4  andMw  <  4,  respectively.  At  these  frequencies,  where  dominant 
wavelengths  are  tens  of  kilometers,  we  assume  a  point  source  for  the  low-magnitude  regional  events  investigated  in 
this  study.  The  point  source  assumption  allows  for  linearization  in  the  time-domain,  which  is  where  we  carry  out  the 
least-squares  inversion. 


RESEARCH  ACCOMPLISHED 


We  analyzed  events  that  were  digitally  recorded  with  a  high  signal-to- 
broadband  stations.  Three-component  data  were  collected  from  a  total 
Network,  IRIS/USGS,  Berkeley  Digital  Seismic  Network,  Trinet,  and 
(LLNL)  network  (Figure  1).  All  data  is  freely  available  from  IRIS  via 
data,  which  is  available  on  compact  disk  (Walter  et  al.,  2004).  Not  all 
1 6  stations  were  used  in  the  inversion  of  the  explosion  data,  which 
are  listed  in  Figure  1.  We  remove  the  instmment  response,  rotate  to 
the  great-circle  frame,  integrate  to  obtain  displacement,  and  use  the 
same  filter  as  for  the  synthetic  seismograms.  The  LLNL  network 
(white  triangles  in  Figure  1)  was  composed  of  Sprengnether 
instmments  with  limited  long-period  response,  and  for  those  data 
we  used  a  passband  of  10-30  seconds  for  both  the  data  and 
synthetics. 

We  calibrated  the  algorithm  by  calculating  the  full  moment  tensor 
for  the  1992  Little  Skull  Mountain  event.  We  find  a  solution  at  all 
depths  within  5  km  of  the  reported  depth.  The  depth  that  produces 
Green’s  functions  that  best  fit  the  data  is  used  in  the  final  solution. 

Fit  is  quantified  by  the  variance  reduction  (VR),  which  is  a 
normalized  variance  given  by 
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where  i  are  the  displacements  at  all  times  for  all  components  at  all 
stations,  and  d  and  ^  are  the  data  and  synthetic,  respectively. 

We  also  allow  the  Green’s  hinctions  calculated  at  a  given  distance 
to  shift  relative  to  the  data  to  address  small  hypocentral  errors  and 
uncertainty  in  the  velocity  model  used  to  compute  the  Green’s 
functions.  The  shift  that  produces  the  best  fit  is  used  in  the  final 
solution.  We  limit  the  shift  to  less  than  5  and  3  sec  for  high-pass 
comers  of  0.05  and  0.10  Hz,  respectively.  The  allowed  time  shift  is 
large  enough  to  make  up  for  small  hypocentral  errors,  but  small 
enough  to  disallow  cycle-skipping  that  could  produce  erroneous 
mechanisms.  The  sensitivity  of  the  time  shift  relative  to  the 
assumed  velocity  model  will  be  discussed  later  in  the  paper.  The 
full  moment  tensor  solution  is  decomposed  to  an  isotropic  and 
deviatoric  component.  We  calculate  the  total  scalar  moment  (Mq)  as 
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Figure  1.  Map  of  the  Western  US  with 

stations  (blue  inverted  triangles), 
earthquakes  (yellow  stars), 
explosions  (red  stars),  and 
collapses  (green  stars)  used  in  this 
study.  The  bottom  panel  is  a  blow¬ 
up  of  the  Nevada  Test  Site  (NTS) 
region  with  the  NTS  outlined  in 
black  and  in  the  top  panel  in  red. 
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a 


b 


Figure  2.  Source-type  plot  for  the  Little  Skull  Mt.  earthquake  (dark  grey  circle),  NTS  test  HOYA  (light  grey 
diamond),  and  bootstrap  population  of  the  NTS  test  DIVIDER  (black  dots)  along  with  its  95% 
confidence  region  (grey  ellipse),  a)  The  source-type  parameters  (A^,  -Is)  given  on  a  linear  plot,  b) 

The  source-type  plot  of  Hudson  et  al.  (1989)  with  theoretical  mechanisms  plotted  as  well. 

defined  by  Bowers  and  Hudson  (1999),  where  Mq  is  equal  to  the  sum  of  the  isotropic  moment  (Mjso  = 
(Mii+M22+M33)/3)  and  deviatoric  moment  (Mdev)?  which  is  the  largest  deviatoric  eigenvalue.  For  the  Little  Skull 
Mountain  event  we  find  Mq  =  3.7  x  10^^  N-m  (Mw  =  5.6),  and  the  solution  has  a  negligible  isotropic  moment  (Mjso  = 
-0.31  X  10^^  N-m)  so  there  is  little  change  between  the  hall  and  deviatoric  solutions.  The  solution  fits  the  data  very 
well  (Figure  2b)  and  is  similar  to  the  double-couple  solution  of  Walter  (1993),  the  deviatoric  solution  of  Ichinose  et 
al.  (2003),  and  the  full  solution  of  Dreger  and  Woods  (2002).  With  the  same  algorithm  we  calculate  the  full  moment 
tensors  of  17  nuclear  test  explosions  at  the  Nevada  Test  Site  (NTS)  (Figure  1).  In  the  case  of  explosions  and 
collapses  we  calculate  Green’s  functions  at  a  depth  of  1  km.  The  sensitivity  of  this  assumption  will  be  investigated 
later  in  the  paper.  An  example  of  the  analysis  is  given  by  the  solution  for  the  1991  HOYA  test,  where  the  largest 
component  in  the  decomposition  is  isotropic,  and  it  contributes  70%  of  the  total  scalar  moment. 

Source-Type  Plot 

It  is  difficult  to  grasp  the  source-type  from  the  standard  focal  mechanism  plot  for  events  with  a  large 
non-double-couple  (DC)  component.  For  example,  one  cannot  discern  the  relative  contributions  of  the  isotropic  and 
deviatoric  components  from  the  full  focal  mechanism  for  the  HOYA  explosion.  In  order  to  get  at  the  tectonic 
contribution  to  the  explosion,  one  could  separate  the  deviatoric  component  into  a  DC  and  a  compensated  linear 
vector  dipole  (CLVD)  that  share  the  orientation  of  the  major  axis,  but  decompositions  of  this  type  are  non-unique, 
where  for  example  the  DC  and  CLVD  decomposition  could  be  replaced  by  two  DCs  (see  Julian  et  al.  [1998])  for 
further  decompositions).  In  an  attempt  to  better  characterize  mechanisms,  we  follow  the  source-type  analysis 

—fYl’  M 

described  in  Hudson  et  al.  (1989)  and  calculate  -2s  and  k,  which  are  given  by  ^  ^  and  k  -  - ^ — r,  where 

Kl  l^iso|+KI 

m'l,  m'2  and  m's  are  the  deviatoric  principal  moments  for  the  N,  P,  and  T  axes,  respectively,  and  \m'i\  <  \m'2\  < 

A  measure  of  the  departure  of  the  deviatoric  component  from  a  pure  DC  mechanism,  s,  is  0  for  a  pure 
double-couple  and  ±0.5  for  a  pure  CLVD.  A  measure  of  the  volume  change  is  k,  where  ±1  would  be  a  full  explosion 
and  -  la  full  implosion.  For  the  Little  Skull  Mountain  earthquake  and  the  NTS  explosion,  HOYA,  -2  s  and  k  are 
given  in  Figure  2a.  The  earthquake  is  almost  at  the  origin,  which  defines  a  pure  DC,  whereas  the  nuclear  test  is  near 
where  a  theoretical  explosion  would  plot.  In  order  to  estimate  formal  error  in  the  fit,  we  create  moment  tensor 
populations  by  bootstrapping  the  residuals  of  the  fit  ^times  with  replacement  and  then  use  those  populations  of  size 
n  to  calculate  -2^  and  k,  resulting  in  their  own  populations  to  which  we  fit  normal  distributions.  Figure  2a  shows  the 
population  of  ^  =  1,000  along  with  the  95%  confidence  region  for  the  DIVIDER  explosion.  Increasing  n  resulted  in 
no  change  to  the  confidence  regions. 

Hudson  et  al.  (1989)  transform  the  parameters  -2s  and  k  so  that  the  displayed  plot  will  have  equal  normal 
probability  areas  based  on  the  assumption  that  the  smallest  principal  moments  can  take  any  value  between  ±  the 
largest  absolute  principal  moment  (Julian  et  al.,  1998).  The  plot  derived  this  way  is  the  source-type  plot,  and  it  is 
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shown  in  Figure  2b  for  the  parameters  from  the  Little  Skull  Mt.  earthquake  and  HOYA  explosion.  Figure  2b  also 
shows  the  transformed  bootstrap  population  for  the  DIVIDER  explosion  and  its  associated  95%  confidence  region. 
The  transformation  makes  the  assumption  of  normality  in  the  error  distribution  valid  as  can  be  seen  by  the  improved 
fit  of  an  error  ellipse  to  the  bootstrap  population  between  Figure  2a  and  b.  The  Hudson  et  al.  (1989)  plot  is  a  superior 
way  to  display  source-type  and  analyze  error  in  the  parameters.  The  error  ellipses  are  not  shown  for  the  Little  Skull 
Mt.  earthquake  or  HOYA  explosion  examples  because  the  error  regions  are  too  small  to  notice  a  difference  due  to 
the  transformation. 

We  carry  out  similar  analyses  for  1 1  more  earthquakes  and  three  collapses  (one  cavity  and  two  mine)  and  produce 
the  source-type  plot  in  Figure  3  along  with  the  95%  confidence  regions.  The  nuclear  tests  occupy  the  region  where 
k>  0.5,  the  earthquakes  cluster  near  the  origin,  and  the  collapses  plot  almost  exactly  at  (l,-5/9),  which  is  the  location 
for  a  closing  crack  in  a  Poisson  solid.  Deviations  from  these  trends  will  be  discussed  later. 

Sensitivity  Analysis 

The  relatively  small  area  of  the  confidence  regions  given  in  Figure  3  and  the  excellent  synthetic  seismogram  fit  to 
the  data  as  quantified  by  VR  gives  us  great  confidence  that  the  assumed  velocity  model  and  depth  are  correct  and  the 
estimated  moment  tensor  solutions  are  robust.  However,  these  measures  of  goodness-of-fit  assume  the  underlying 
model  used  to  invert  the  data  is  correct.  In  the  following  section  we  will  test  these  assumptions  with  synthetic  data 
from  a  theoretical  explosion  (-26=0,  k^\)  created  for  two  experimental  geometries.  The  first  geometry,  referred  to 
as  ‘Ideal,’  is  eight  stations  at  distance  increments  between  100  and  300  km  each  separated  by  45°.  The  second 
station  geometry  mirrors  the  analysis  for  the  HOYA  explosion.  The  synthetic  data  are  filtered  in  the  same  two  bands 
(20-50s  and  10-50s)  used  in  the  analysis  and  when  combined  with  the  two  geometries  results  in  4  scenarios. 


Explosion 


Figure  3.  Source-type  plot  of  the  12  earthquakes  (blue),  17  explosions  (red),  3  collapses  (green),  and  their 

associated  95%  confidence  regions  (shaded).  The  magnitude  of  the  event  is  given  by  the  symbol.  The 
abscissa  measures  the  amount  of  volume  change  for  the  source  and  the  ordinate  measures  the 
departure  from  pure  DC.  Theoretical  mechanisms  (crosses)  are  plotted  for  comparison. 
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The  error  analysis  presented  above  is  due  to  a  misfit  of  the  data  by  the  least-squares  inversion.  Part  of  the  misfit  may 
be  due  to  nonstationary  noise,  and  we  test  the  sensitivity  of  the  inversion  to  different  SNRs.  In  order  to  best 
approximate  real-world  noise  conditions,  we  derive  the  noise  signal  from  data  prior  to  the  first  arrival  from  the 
nuclear  test  METROPOLIS  (10  March  1990)  at  station  ANMO  for  all  three  components.  The  amplitude  of  this  noise 
signal  is  bandpassed  to  match  the  synthetic  data  and  multiplied  by  a  factor  so  as  to  create  a  final  synthetic  signal 
with  the  desired  SNR  (ratio  of  synthetic  data  root-mean-square  amplitude  to  noise  root-mean-square  amplitude). 


The  noise  analysis  has  very  little  frequency  dependence,  so  for  clarity  we  only  show  results  from  the  analysis  in  the 
20-50  sec  frequency  band  in  Figure  4a.  The  Ideal  configuration  produces  the  best  scenario  where  a  large  k  is 
retrieved  (>0.3)  when  the  SNR  is  greater  than  2.  For  all  scenarios  k  >  0.5  when  SNR  >  5.  Typically,  we  use  data 
with  an  SNR  greater  than  10;  however,  there  are  a  few  cases  where  the  SNR  is  close  to  3.  An  example  of  this  type  of 
data  is  for  the  DIVIDER  explosion,  which  produced  a  signal  that  was  right  on  the  limit  of  acceptable  SNR  but  still 
produced  a  well-fit  solution. 


Another  source  of  error  not  incorporated  into  the  formal  error 
analysis  is  incorrectly  calculated  Green's  functions  due  to 
ignorance  of  the  true  event  depth.  The  method  that  produces  the 
results  presented  above  attempts  to  find  an  optimal  depth  for  the 
earthquakes  by  perturbing  the  reported  depth  a  few  kilometers, 
performing  the  inversion,  and  finding  the  best-fit  solution.  For  all 
explosions  and  collapses,  the  depth  is  fixed  at  1  km.  If  this 
method  were  to  be  used  for  an  event  with  an  unknown  source 
type,  the  depth  could  be  an  important  source  of  error,  as  well  as  an 
important  parameter  for  identification.  We  perform  another 
synthetic  test  in  which  an  explosion  at  1  km  is  inverted  with 
Green's  functions  calculated  at  varying  depths. 

The  source  depth  analysis  is  not  greatly  affected  by  the  two 
station  configurations  considered  here,  therefore  we  only  show 
results  for  the  HOYA  configuration  in  Figure  4b.  The  result  at  an 
incorrect  depth  of  2  km  is  virtually  indistinguishable  from  the  true 
answer.  When  the  source  is  moved  to  3  km  depth  there  is  a  small 
step  decrease  in  k  due  to  a  layer  in  the  velocity  model  that  begins 
at  a  2.5-km  depth.  However,  k  >  0.5  for  incorrect  depths  <17  km 
with  slightly  more  sensitivity  in  k  and  a  worse  fit  in  the 
high-frequency  band  (10-50  sec)  compared  to  the  low-frequency 
band  (20-50  sec).  The  relative  insensitivity  of  the  solution  to 
mislocated  depth  for  an  explosion  is  different  than  is  observed  for 
DC  events.  Dreger  and  Woods  (2002)  show  that  the  VR  of  the 
Little  Skull  Mountain  earthquake  solution  is  definitively 
maximized  at  the  assumed  true  event  depth.  Thus,  while  the  depth 
sensitivity  of  explosions  is  poor,  the  method  is  able  to  determine 
the  depth  of  non-explosion  sources,  which  also  provides  an 
important  level  of  event  screening. 

Finally,  we  test  how  error  in  the  assumed  Earth  structure  is 
mapped  through  the  Green's  functions  to  error  in  the  solution.  We 
start  with  the  well-calibrated  Song  et  al.  (1996)  velocity  model 
(Table  1)  and  perturb  the  velocities  and  depths  of  the  layers  using 
averaged  parameters  from  another  plausible  velocity  model 
(WestUS;  Ammon,  1999)  and  a  model  from  Southern  California 
(SoCal;  Dreger  and  Helmberger,  1990). 

In  order  to  produce  a  sensitivity  test  that  best  mimics  our  analysis, 
we  use  the  time  shift  rule  to  filter  the  models.  This  means  that  we 
only  allow  velocity  models  that  produce  Green’s  functions  where 
the  time  shift  between  data  and  synthetics  that  produces  the  best- 
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Figure  4.  Sensitivity  analysis,  a)  Noise  is 

added  to  the  inversion  of  20-50  sec 
synthetic  data  while  velocity  model 
and  depth  (1  km)  are  kept  fixed  for 
the  HOYA  (circle)  and  Ideal 
(triangle)  scenarios,  b)  The  inversion 
using  the  HOYA  configuration  is 
carried  out  assuming  an  incorrect 
depth  while  the  velocity  model  is 
kept  fixed  for  data  in  the  20-50-sec 
(circle)  and  10-50-sec  (triangle) 
band,  c)  The  inversion  using  the 
HOYA  configuration  for  20-50  sec 
synthetic  data  is  carried  out  for 
different  three-layer  velocity  models 
where  the  data  are  not  shifted 
relative  to  the  Green’s  functions  (left 
panel,  circles)  and  allowed  to  shift 
less  than  5  sec  (right  panel, 
triangles).  The  symbols  are  colored 
as  a  function  of  VR. 
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Figure  5.  Vanishing  traction  sensitivity.  Synthetic  data  for  a  pure  explosion  (A^=l)  is  inverted  at  depths  less 

than  1  km  for  varying  SNR  and  the  four  scenarios  discussed  in  the  text,  a-d)  Resolved  Miso  for  SNR 
values  of  2  (circle)  6  (inverted  triangle)  and  10  (triangle)  where  the  value  for  an  inversion  without 
noise  (SNR=oo)  is  given  by  the  black  line  (100%).  k  is  given  by  the  color,  e-h)  Resolved  Mdev  for  SNR 
values  of  2,  6,  and  10  with  the  same  symbols,  where  the  total  scalar  moment  for  an  inversion  without 
noise  (SNR=oo)  is  given  by  the  black  line  (100%),  and  Mdev  should  be  0.  -28  is  given  by  the  color. 


fit  solution  is  less  than  or  equal  to  5  or  3  sec  from  the  theoretical  arrival  time  for  high-pass  comers  of  0.05  or  0.10 
Hz,  respectively.  Primarily  due  to  the  velocity  model  filtering  there  is  little  difference  among  the  scenarios  so  we 
only  show  source-type  plots  for  the  HOYA  configuration  in  the 

20-50-sec  frequency  band  in  Figure  4c.  For  this  scenario  the  number  of  acceptable  models  is  reduced  to  88,  and 
although  not  all  possible  combinations  of  model  parameters  are  used,  each  parameter  perturbation  given  is  employed 
at  least  once. 


Without  shifting,  there  are  a  few  velocity  models  that  produce  well-fit  solutions  (VR>90%)  with  mechanisms  that 
are  almost  purely  DC.  However,  when  shifting  is  allowed,  all  velocity  models  produce  good  fits  with  highly 
explosive  sources  (A:~>0.4). 


Free-Surface  Effects 


Figure  6.  Map  of  the  YSKP  with  the  North 
Korea  test  (red  star)  and  nearby 
earthquake  (yellow  star)  as  well  as 
the  stations  used  in  this  study  (blue 
triangles). 


Another  consideration  is  the  ability  to  resolve  displacements  for 
explosions  near  the  surface.  Since  tractions  normal  to  the  vertical 
vanish  at  the  free  surface,  the  excitation  coefficients  associated 
with  those  tractions  must  vanish  (Julian  et  al.,  1998).  Therefore,  at 
the  free  surface  the  moments  of  Mi  3,  M23,  and  the  isotropic  part  of 
the  Mij  cannot  be  resolved.  Given  and  Mellman  (1986)  showed 
that  at  a  source  depth  of  1  km,  the  fundamental  mode  excitation 
functions  associated  with  the  moments  listed  previously 
effectively  go  to  zero.  We  investigate  the  potential  problems 
associated  with  vanishing  traction  at  the  free  surface  by  inverting 
noisy  data  from  a  synthetic  explosion  source  at  depths  between 
200  and  1,000  m  in  a  three-layer  ID  velocity  model  using  Green’s 
functions  calculated  at  those  same  depths. 

The  ability  to  resolve  an  explosive  component  is  dependent  on  the 
station  distribution,  frequency,  and  SNR  of  the  analysis,  therefore 
Figure  5  shows  all  four  scenarios.  An  explosive  component 
(A:>0.5)  can  be  resolved  under  favorable  noise  conditions  at  a 
depth  greater  than  300  m  for  all  scenarios,  though  with  error  in 
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•  99%VRfit 

Figure  7.  Network  sensitivity  solutions  for  the  North  Korea  explosion  (left)  and  earthquake  (right).  The 

actual  solution  is  given  by  the  stars  and  other  solutions  for  a  given  threshold  VR  are  shown.  The 
95%  VR  fit  is  given  with  (dark  gray)  and  without  station  BJT  (light  gray). 

Miso  between  50-150%  (Figure  5a-d).  The  error  is  inversely  proportional  to  the  depth.  For  all  scenarios,  but  the 
HOYA  configuration  at  20-50  sec  (Figure  5a),  favorable  noise  means  SNR  >  6.  The  change  in  Mjso  is  due  to  a 
change  in  M33  relative  to  the  other  dipole  components,  and  this  produces  an  erroneous  deviatoric  component.  The 
moment  of  deviatoric  component  can  be  up  to  50%  of  the  theoretical  isotropic  moment  (Figure  5e-h)  and  since  it  is 
related  to  the  error  in  Mjso  it  is  inversely  proportional  to  the  depth.  At  less  than  200  m  depth,  the  synthetic 
displacements  become  too  small  and  the  solution  using  these  particular  Green’s  functions  is  unreliable. 

Network  Sensitivity  Solution 

We  extend  the  database  of  solutions  by  incorporating  the  North  Korea  test  and  a  nearby  earthquake  that  occurred  on 
16  December  2004  (Figure  6).  We  employ  the  records  of  four  stations  that  recorded  the  events  well  in  the  period 
band  of  interest.  We  use  the  solutions  to  present  a  new  way  to  estimate  constraints  on  source-type.  As  discussed 
above,  the  ability  to  resolve  a  well-constrained  solution  is  dependent  on  station  configuration,  data  bandwidth,  and 
SNR.  The  analysis  presented  in  the  previous  section  tried  to  assess  error  in  those  experimental  conditions,  but  it  is 

+v  +v 


O  Example  shown 


Figure  8.  Network  sensitivity  solutions  for  the  North  Korea  explosion  station  configuration  and  Green’s 

functions  (left)  and  earthquake  setup  (right).  The  best  solution  is  given  by  the  stars  and  different 
levels  of  fit,  as  quantified  by  VR,  are  shown.  The  source-types  of  the  solutions  given  in  Figures  10 
and  11  are  also  shown  (blue  circles). 
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difficult  to  state  steadfast  rules  for  what  source-types  can  be  resolved  for  all  conditions,  when  different  conditions 
lead  to  different  levels  of  confidence  in  the  solution.  Therefore,  we  calculate  two  types  of  forward  confidence 
analyses,  which  we  call  network  sensitivity  solutions. 

The  first  type  of  solution  assumes  a  very  large  SNR  and  produces  explosive  or  double-couple  synthetics  for  a  given 
station  configuration  using  the  Green’s  function  employed  in  the  inverse  solution.  In  this  case,  the  MDJ2  model  is 
used  to  calculate  the  Green’s  functions.  These  model  explosions  and  earthquakes  are  then  compared  with  synthetic 
data  produced  by  thousands  of  random  seismic  moment  tensors  that  represent  a  uniform  distribution  of  all  possible 
moment  tensors.  A  VR  is  calculated  between  the  synthetics  produced  by  the  distribution  of  moment  tensor  and  the 
perfect  explosion  or  DC.  In  this  way,  we  produce  a  source-type  plot  that  shows  solutions  for  a  given  threshold  VR. 
Figure  7  shows  the  source-type  plot  constructed  in  this  way  for  the  North  Korea  test  and  nearby  earthquake 
configurations.  Without  station  BJT,  most  solutions  fit  the  model  explosive  just  as  well  as  an  explosion 
(Figure  7,  left),  whereas  any  source  with  an  isotropic  source  greater  than  50%  does  not  fit  the  DC  synthetics  well 
(Figure  7,  right).  The  use  of  this  type  of  network  sensitivity  solution  can  give  the  user  a  feeling  for  what  level  of 
confidence  can  be  expected  in  a  given  solution  for  the  experimental  setup  employed. 

The  network  sensitivity  solution  is  also  produced  using  the  actual  data  in  place  of  the  synthetic  explosion  or  DC. 
This  plot  can  help  the  user  understand  the  limits  of  confidence  in  the  best  solution.  Figure  8  shows  this  type  of 
network  sensitivity  solution  for  the  same  events  as  in  Figure  7.  The  solutions  are  fairly  well-constrained,  but  there 


Strike  =  42/305  40/305 
Rake  =  -150/- 14 -149/- 1  1 
Dip  =  78/61  80/60 
M0=  l.le-i-22  1.2e-h22 
MW  =  4.0  4.0 
DC  (%)  =  88  SO 
CLVD(%)  =  6  1 
ISO(%)  =  6  19 
VR(%)  =  73.7  70.2 


Figure  9.  Waveform  fits  and  focal  mechanism  for  the  best-fit  (red)  and  fairly-well  fit  (blue)  solutions  using 
data  from  the  earthquake  near  the  North  Korea  test.  The  station  name,  azimuth,  distance,  and 
maximum  displacement  in  cm  are  given  next  to  the  traces.  Information  on  the  sources  is  also  given. 
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are  some  possible  solutions  with  good  VR  that  are  near  the  DC  section  of  the  source-type  plot,  for  the  explosion  data 
(Figure  8,  left),  and  solutions  that  have  some  isotropic  component,  for  the  earthquake  data  (Figure  8,  right).  An 
understanding  of  what  the  VR  means  for  these  potentially  confidence-busting  solutions  can  be  acquired  by 
examining  the  data- fit.  Figures  9  and  10  give  the  data  fit  for  the  best  solutions  (stars  in  Figure  8)  and  those  for 
‘problematic’  solutions  with  a  high  VR  (blue  circles  in  Figure  8).  Visual  inspection  of  the  fits  show  that  they  are 
virtually  indistinguishable  and  that  both  solutions  may  be  possible  given  the  data  limitations.  This  is  not  so  much  a 
problem  for  the  earthquake  solution  (Figure  9)  since  the  percent  DC  is  still  relatively  high,  but  the  nuclear  test 
solutions  are  fairly  different  (Figure  10). 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  populations  of  earthquakes,  explosions,  and  collapses  separate  in  the  source-type  plot.  These  initial  results  are 
very  encouraging  and  suggest  a  discriminant  that  employs  the  source-type  plot  parameters  (-2^,  k).  Another 
advantage  of  the  source-type  plot  is  its  display  of  2-D  error  regions.  In  this  way  one  can  test  a  hypothesis  that  an 
event  has  a  non-DC  component.  The  source-type  analysis  can  also  be  utilized  to  estimate  model-based  error  as  well. 
The  error  introduced  by  ignorance  of  the  event  location  and  Earth  structure  can  be  calculated  with  a  Monte  Carlo 
approach,  where  several  solutions  are  computed  for  a-priori  distributions  of  the  hypocentral  location  and  Earth 
model  obtained  from  independent  analyses. 

We  try  to  give  some  insight  to  the  depth  sensitivity  of  the  method  with  Figure  4b,  which  shows  that  the  use  of  this 
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Figure  10.  Waveform  fits  and  focal  mechanism  for  the  best-fit  (red)  and  fairly-well  fit  (blue)  solutions  using 
data  from  the  North  Korea  test.  See  Figure  9  for  more  information. 
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method  as  a  precise  depth  discriminant  is  not  plausible  for  the  frequencies  used  here,  though  sensitivity  does 
increase  for  the  higher-frequency  band.  These  results  are  a  demonstration  of  the  fact  that  an  isotropic  radiation 
pattern  has  no  sensitivity  to  takeoff  angle,  which  depends  on  depth.  As  shown  by  Dreger  and  Woods  (2002),  there  is 
limited  resolution  of  the  shallow  depth  of  explosions  using  regional  distance  data.  Although  an  explosive  radiation 
pattern  alone  does  not  have  depth  sensitivity,  the  relative  excitation  of  low-frequency  body  waves  (Pnl)  and 
Rayleigh  waves  does  enable  the  method  to  discern  the  relatively  shallower  depths  of  explosions  compared  to 
earthquakes. 

The  velocity  model  analysis  shown  in  Figure  4c  suggests  that  the  maximum- shift  rule  used  in  the  analysis  is  a  good 
proxy  for  evaluating  the  appropriateness  of  the  velocity  model.  The  level  of  departure  of  a  given  velocity  model 
from  the  true  model  is  station  distribution,  frequency,  and  SNR  dependent.  Therefore,  it  is  a  good  idea  to  perform 
this  style  of  sensitivity  test  to  evaluate  the  amount  of  deviation  a  certain  experimental  setup  will  allow,  because  if 
the  velocity  model  is  poorly  calibrated,  then  a  good  fit  to  the  data  can  be  obtained,  but  the  solution  may  be 
inaccurate.  The  network  sensitivity  solution  is  a  good  way  to  accomplish  this  type  of  testing. 

The  change  in  moment  due  to  the  loss  of  traction  at  the  free  surface  affects  yield  estimation,  though  event 
discrimination  is  still  reliable  at  high  SNR.  A  result  of  this  change  in  moment  is  that  the  deviatoric  moment  becomes 
non-zero  and  could  be  significant  at  very  shallow  depths  (Z<500  m)  and  low  SNR  (SNR<6).  The  moment  manifests 
as  a  CLVD  component,  which  means  that  interpretation  of  non-isotropic  energy  may  be  flawed  for  shallow  events 
even  with  high  SNR  data.  Though,  as  Figure  5  suggests,  this  effect  is  station  configuration,  frequency,  and  SNR 
dependent. 
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